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Abstract 

Chiral symmetry breaking is calculated as a function of cooling in quenched 
lattice QCD. A non-zero signal is found for the chiral condensate beyond 
one hundred cooling steps, suggesting that there is chiral symmetry breaking 
associated with instantons. Quantitatively, the chiral condensate in cooled 
gauge field configurations is small compared to the value without cooling. 
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Typeset using REVTgX 



Recently Chu et al [|TJ presented evidence from a quenched SU(3) lattice gauge theory that 
instantons play a dominant role in determining hadronic correlation functions. Motivated 
by this work, a calculation was done in quenched SU(2) lattice gauge theory exploring the 
effect of cooling ||. In Ref. || it was found that the spin-dependent potential and the chiral 
symmetry breaking order parameter disappear very quickly upon cooling. In a subsequent 
nonquenched SU(3) calculation by Edwards et al it was found that cooling has a large 
effect on the hadron spectrum, tending, for example, to diminish spin-slittings. This supports 
the conclusions of Ref. |§. However, Edwards et al |§ did not see the rapid disappearance of 
the chiral order parameter found in the quenched SU(2) calculation. For this reason we have 
decided to revisit the problem of chiral symmetry breaking and cooling. A calculation was 
done in quenched SU(3) lattice gauge theory to see if the different chiral symmetry breaking 
pattern observed in Ref. and in Edwards et al is due to the use of SU(2) versus SU(3) 
color or due to the quenched versus nonquenched nature of the simulation. 

In this paper we present results for the chiral order parameter calculated in a set of 
quenched SU(3)-color gauge configurations to which local cooling ||,[|] has been applied. A 
persistent non- vanishing signal for chiral symmetry breaking was found even for cooling of a 
hundred sweeps or more. This differs from what was observed in the earlier quenched SU(2) 
simulation but is consistent with the finding of Edwards et al || in their nonquenched 
SU(3) calculation. Quantitatively, after cooling, the chiral order parameter is small (in 
lattice units) compared to the uncooled value. 

The calculations were done using standard methods of Ref. carried over to SU(3). 
The fermion action is 

s f = \ u ^ x ) x( x + A) - xi x + fr) u t( x ) x(x)] + Yl m x( x )x( x ), 

X,fl X 

= XM({U}) X , (1) 

where x, % are single- component fermion fields, rj^x) is the staggered fermion phase ||, m 
is the mass (in units of the inverse of the lattice spacing a) and the U's are gauge field links. 
Antiperiodic boundary conditions were used for the staggered fermion fields in all directions. 
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The chiral symmetry order parameter is calculated from the inverse of the fermion matrix 
M of Eq. (g) 

(XX) = ^(^M- 1 ({U})}, (2) 

where V is the lattice volume and the angle brackets indicate a gauge field configuration 
average. The random source method Jrj was used to calculate Tr.M _1 ({i7}). Twelve 
Gaussian random sources were used for each gauge field configuration. 

Calculations were done on a 10 4 lattice at (3 = 5.7 and on a 12 at (3 — 5.9. The SU(3) 
pseudo-heatbath method was used for updating the quenched gauge field configurations. 
After 4000 sweeps to thermalize, configurations were analyzed every 300 sweeps. A total of 
30 configurations at f3 — 5.7 and 20 configurations at j3 = 5.9 were included in the analysis. 
Each configuration was cooled using a local cooling algorithm [|]|| and (xx) was calculated 
along the way. 

The results are tabulated in Table |. The results for (xx) as a function of quark mass 
at no cooling and after 20 cooling sweeps are shown in Fig. [I] and Fig. ^| for (3 = 5.7 and 
(3 = 5.9 respectively. The value of (xx) obtained in a free field theory is also plotted in the 
figures. 

Following Bowler et al |§, the values at zero quark mass (extrapolated values) were 
obtained by fitting to the function 

(xx)(m) = (xx)o + (xx)i m + (xx)a m 2 ■ (3) 

The resulting chiral condensate values (xx)o are plotted in Fig. versus the number of 
cooling sweeps. A clear nonzero signal for chiral symmetry breaking is obtained confirming 
the conclusion of Edwards et al . This is in clear contrast to the results of Ref . (see Figs. 
9-11) where the chiral order parameter was found to disappear very quickly with cooling. It 
is fair to say, however, that the value of the chiral order parameter found in our calculation, 
although nonzero after cooling, is small compared (in lattice units) to the uncooled value. 
The behavior of topological properties under cooling is well known |§. Here we 
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give a brief description to put our calculation into context. The simple transcription of 
F^ u {x)F^ u {x) to the lattice |L(J is used for the topological charge density 



where U^, U pa are plaquettes of gauge field links. Fig. |] shows the expectation values of 
topological charge squared Q 2 (Q = J2x1t( x )) anc ^ a measure of the topological activity, 
At = \qt(x)\. Note that without cooling the topological activity A t is very large due to 
short distance fluctuations in qt(x). The numbers are too large to be plotted in Fig. f|(b) and 
are 45.02 ± 0.05 at (3 = 5.7 (10 4 lattice) and 87.6 ± 0.1 at /3 = 5.9 (12 4 lattice). In the first 
ten or so cooling steps, short distance fluctuations disappear very quickly, the topological 
activity drops rapidly and (Q 2 ) takes a value which changes only slowly with more cooling. 
In the region up to about 60 cooling steps {A t ) 2 is larger than (Q 2 ), indicating that one 
does not yet have gauge field configurations consisting of single isolated instantons. Rather, 
the vacuum in this region of cooling is sometimes characterized in terms of "overlapping" 
or "interacting" instantons |1]||. Beyond about 60 cooling steps (A t ) 2 < (Q 2 ) as one would 
expect for a sample of instanton configurations with differing topological charge. The rate 
at which the topological activity changes becomes very small as instantons, which are not 
absolutely stable on the lattice, slowly shrink and fade away. By 100 cooling steps one 
expects to be in the domain of a vacuum dominated by isolated instantons [0,0,0. Our 
finding that the value of the chiral order parameter is non-zero and essentially constant 
beyond 60 cooling steps, we would say, gives an indication that, in an SU(3) gauge theory, 
instantons produce chiral symmetry breaking, although it would not suggest that instantons 
dominate chiral symmetry breaking. 

At present we have no explanation for the different behavior observed in the SU(2) 
and SU(3) calculations. One possibility is that the SU(2) calculation, which was done at 
P = 2.4, is just deeper in the weak coupling region than the SU(3) calculations done here at 
(3 = 5.7 and 5.9. Indeed, comparing the SU(3) calculations at the two f3 values may suggest 
a diminished persistence of chiral symmetry breaking at (3 — 5.9 but this conclusion can not 




(4) 



4 



be made definitively due to the smallness of the chiral condensate (and concomitant large 
uncertainty) at the larger (3 value. It may also be that there is still some more fundamental 
difference between the SU(2) and SU(3) theories which has yet to be found. 

It is a pleasure to thank H.D. Trottier for helpful discussions. This work was supported 
in part by the Natural Sciences and Engineering Research Council of Canada. 
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TABLES 

TABLE I. Chiral order parameter values calculated for different values of the quark mass along 
with the value extrapolated to zero mass. 



ma 





Cooling step 


0.2 


0.15 


0.1 


0.05 


Extrapolated value 


5.7 





0.889(2) 


0.788(2) 


0.660(2) 


0.494(3) 


0.30(1) 




20 


0.3659(3) 


0.2835(5) 


0.1986(7) 


0.115(2) 


0.028(4) 




40 


0.3613(5) 


0.2783(5) 


0.1930(8) 


0.110(2) 


0.023(5) 




60 


0.3600(4) 


0.2768(6) 


0.1913(9) 


0.109(2) 


0.020(4) 




80 


0.3593(5) 


0.2761(6) 


0.1907(9) 


0.107(2) 


0.020(4) 




100 


0.3592(5) 


0.2757(6) 


0.1899(9) 


0.107(2) 


0.020(5) 


5.9 





0.7775(7) 


0.6589(8) 


0.5104(9) 


0.316(1) 


0.091(4) 




20 


0.3595(2) 


0.2760(2) 


0.1892(3) 


0.1009(5) 


0.009(2) 




40 


0.3583(2) 


0.2745(2) 


0.1874(3) 


0.099(6) 


0.007(2) 




60 


0.3576(2) 


0.2738(2) 


0.1866(3) 


0.098(6) 


0.005(2) 




80 


0.3576(2) 


0.2736(2) 


0.1862(3) 


0.098(6) 


0.005(2) 




100 


0.3573(2) 


0.2732(2) 


0.1859(3) 


0.097(6) 


0.005(2) 
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FIG. 1. Plot of the chiral order parameter (in lattice units) as a function of quark mass at 
= 5.7 for no cooling (A) and 20 cooling steps (•). The lattice free-field values are also shown 
). The lines are the extrapolations to zero quark mass using Eq. (||). 
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FIG. 2. Same as in Fig. [I], except for = 5.9. 
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FIG. 3. The chiral order parameter at zero quark mass (xx)o versus number of cooling sweeps 
at (a) = 5.7, and (b) (3 = 5.9. 
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FIG. 4. Expectation value of topological charge squared (a) and topological activity (b) as a 
function of cooling at j3 = 5.7 (■) and (3 = 5.9 (•). 
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